Large collection of astrophysical S-factors and its compact representation 
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Numerous nuclear reactions in the crust of accreting neutron stars are strongly affected by dense 
plasma environment. Simulations of superbursts, deep crustal heating and other nuclear burning 
phenomena in neutron stars require astrophysical S-factors for these reactions (as a function of 
center-of-mass energy E of colliding nuclei). A large database of S-factors is created for about 5,000 
non-resonant fusion reactions involving stable and unstable isotopes of Be, B, C, N, O, F, Ne, Na, 
Mg, and Si. It extends the previous database of about 1,000 reactions involving isotopes of C, O, Ne, 
and Mg. The calculations are performed using the Sao Paulo potential and the barrier penetration 
formalism. All calculated S-data are parameterized by an analytic model for S{E) proposed before 
[Phys. Rev. C 82, 044609 (2010)] and further elaborated here. For a given reaction, the present 
S(-E)-model contains three parameters. These parameters are easily interpolated along reactions 
involving isotopes of the same elements with only seven input parameters, giving an ultracompact, 
accurate, simple, and uniform database. The S{E) approximation can also be used to estimate 
theoretical uncertainties of S{E) and nuclear reaction rates in dense matter, as illustrated for the 
case of the '^*Ne-|-^''Ne reaction in the inner crust of an accreting neutron star. 

PACS numbers: 25.70.Jj;26.50.+x;26.60.Gj,26.30.-k 
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I. INTRODUCTION 

Nuclear burning is an important ingredient of stellar structure and evolution [H-Q . Simulating the evolution and 
observational manifestations of different stars (from main-sequence, to giants and supergiants, presupernovae, white 
dwarfs and neutron stars), requires the rates of many reactions, involving different nuclei - light and heavy, stable 
and neutron-rich. These rates are derived from respective reaction cross sections cf{E), where E is the center-of-mass 
energy of reactants. 

The rates of many reactions which occur in the classical thermonuclear regime - for instance, in main-sequence 
stars - are well known. However nuclear burning in dense plasma of white dwarfs and neutron stars [4] , which affects 
the evolution and many observational manifestations of these objects, may proceed in other regimes, under strong 
effects of plasma screening and pycnonuclear tunneling through Coulomb barrier (e.g., Refs. |5|-|7[). The burning 
powers nuclear explosions in surface layers of accreting white dwarfs (nova events), in cores of massive accreting white 
dwarfs or in binary white dwarf mergers (type la supernovae) [sl-fioj , and in surface layers of accreting neutron stars 
(type I X-ray bursts and superbursts; e.g., Refs. [11-15]). Nova events and type I X-ray bursts are mostly driven by 
the proton capture reactions of the hot CNO cycles and by the rp-process. This burning is thermonuclear, without 
any strong effects of dense plasma environment. Type la supernovae and superbursts are driven by the burning of 
carbon, oxygen, and heavier elements at high densities, where the plasma screening effects can be substantial. It is 
likely that pycnonuclear burning of neutron-rich nuclei (e.g., ^^Ne-|-^*Ne) in the inner crust of accreting neutron stars 
in X-ray transients (in binaries with low-mass companions; e.g. Refs. [l4. ITgI [iTj ) provides an internal heat source for 
these stars. If so, it powers |18 | thermal surface X-ray emission of neutron stars observed in quiescent states of X-ray 
transients (see, e.g., Refs. jl4l. Il9l. [20|) although other energy sources can also be important there (e.g., Ref. [2]|). 

Therefore, in dense plasma of evolved stars, standard classical thermonuclear reaction rates may be unavailable 
and/or inapplicable. Accordingly, one needs to construct many reaction rates, valid over all burning regimes, starting 
from reaction cross sections <j{E). Here we focus on non-resonant fusion reactions involving {Ai,Zi) and (^2,^2) 
reactants {Ai and Zi stand for their mass and charge numbers), which may occur in evolved stars. The cross section 
a{E) is conveniently expressed through the astrophysical S'-factor, 

a{E) = E-^ cxp(-27r77) S{E), (1) 

where = a/{hv) — ^JEr/ E is the Sommerfeld parameter, v = yj2E/ jj. is the relative velocity of the reactants at 
large separations, a = ^1^26^, Er = a^/i/(2?i^) is analogous to the Rydberg energy in atomic physics, and is the 
reduced mass. The factor exp(— 27r7]) is proportional to the probability of penetration through the pure Coulomb 
barrier U{r) — a/r with zero angular orbital momentum, assuming that this pure Coulomb barrier extends to r — >■ 
(as for point-like nuclei); E~^ factorizes out the well-known pre-exponential low-energy dependence of u{E). The 
advantage of this representation is that S{E) is a more slowly varying function of E than (t{E). 

For astrophysical applications, one needs S{E) at low energies, E < a. few MeV. Even for beta-stable nuclei 
experimental measurements of cr{E) at such energies are mainly not available because the Coulomb barrier becomes 
extremely thick making a{E) exponentially small. It will be even more difficult to get such data for neutron-rich 
nuclei. Moreover, the modelling of the processes at different astrophysical sites (such as dynamic reaction network 
modeling of the neutron star crust composition at high densities [22]) requires the knowledge of S'- factors for a large 
variety of reactions, many of which involve neutron-rich nuclei. The experimental study of so many reactions is 
definitely beyond existing and near-future capabilities. Therefore, one must rely on theoretical calculations. 

Previously [23l| we calculated S{E) for 946 reactions involving isotopes of C, O, Ne, and Mg. Also, we proposed 
[23 | a simple analytic model for S{E) and used it to fit calculated S'-factors. The present paper significantly enlarges 
the database. We have extended the calculations to incorporate new reactions between even-even and odd(Z)-even 
isotopes. In particular, we have calculated S{E) between isotopes of Be, B, C, N, O, F, Ne, Na, Mg, and Si. Our 
new database of S{E) contains 4,851 reactions. Moreover, we elaborate and simplify the analytic S'(i?)-model, and 
use it to parameterize all iS-factors with the minimum number of fit parameters producing thus an ultracompact and 
uniform database convenient in applications. 



II. CALCULATIONS 

Consider a set of S'-factors for non-resonant fusion reactions involving various isotopes of 10 elements, Be, B, C, 
N, O, F, Ne, Na, Mg, and Si. The S'-factors for the reactions involving even-even isotopes of C, O, Ne, and Mg have 
been calculated recently in our previous paper (23j . Calculations for other reactions are original and include even-even 
and odd(Z)-even stable, proton-rich, neutron-rich, and very neutron-rich isotopes. Such isotopes can appear during 
nuclear burning in stellar matter, particularly, in the cores of white dwarfs and envelopes of neutron stars. The 
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calculations have been performed using the Sao Paulo potential in the frame of the barrier penetration model |25| . 
Nuclear densities of reactants have been obtained in Relativitic Hartree-Bogoliubov (RHB) theory [2^ employing the 
NL3 parametrization for relativistic mean field Lagrangian [23| and Gogny DIS force for pairing. The model is based 
on the standard partial wave decomposition = 0, 1, . . .) and considers the motion of reacting nuclei in the effective 
potential discussed in Sec. IIII CI [see Eq. (|TOl) there]. The numerical scheme is parameter- free and relatively simple 
for generating a set of data for many non-resonant reactions involving different isotopes. 

The reactions in question are listed in Table HI All reactants are either even-even or odd-even nuclei. We consider 
55 reaction types, such as Si-|-Si and Be-|-B, with the range of mass numbers for both species given in the columns 
2 and 3 of Table H] For each reaction, S{E) has been computed on a dense grid of E (with the energy step of 0.1 
MeV) from 2 MeV to a maximum value -Emax (also given in Table |T|) covering wide energy ranges below and above 
the Coulomb barrier. The last column in Table U presents the number of considered reactions. 

5'(£^)-factors calculated using the Sao Paulo potential have been compared previously [IjH^Hl] with experimental 
data and with theoretical calculations performed using other models such as coupled-channels and fermionic molecular 
dynamics ones. Let us stress that the calculated values of S{E) are uncertain due to nuclear physics effects - due to 
using the Sao Paulo model with the NL3 nucleon density distribution. As shown, for instance, in Ref. [28], typical 
expected uncertainties for the reactions involving stable nuclides are within a factor of 2, with maximum up to a factor 
of 4. For the reactions involving unstable nuclei, typical uncertainties can be as large as one order of magnitude, 
reaching two orders of magnitude at low energies for the reactions with very neutron-rich isotopes. These uncertainties 
reflect the current state of the art in our knowledge of S{E). 



III. ANALYTIC MODEL 



A. General remarks 

Let us recall the main features of the analytic model for S{E) proposed in Ref. and used there to fit the initial 
set of 946 S'-factors. We will elaborate and simplify the model, and fit much larger set of data. 

The model [24] is based on semi-classical consideration of quantum tunneling through an effective potential U{r) 
(which is purely Coulombic at large separations r but is truncated by nuclear interactions at small r when colliding 
nuclei merge). The reaction cross section sd, E < Eq (where Ec is the barrier height) is taken in the form 



<^iE)^— exp 



dr ^2n{U -E) 



(2) 



where ri and r2 are classical turning points; is a slowly varying function of E treated as a constant; it has the 
same dimension as S{E), but should not be confused with it. 

To obtain tractable formulae for S{E), in Ref. [24| we employed the natural and simplest approximation of U{r), 



U{r) 



a 



at r > i?ci. 



U{r) = Ec 



1-/3 



[r - Ref 



at r < Rci, 



(3) 



which is a pure Coulomb potential at r > Rci and an inverse parabolic potential at smaller r (see Fig. 1 in Ref. [2J|). 
The parabolic segment truncates the effective interaction at small separations; Ec — U{Rc) is the maximum of U{r), 
that is the barrier height. We required U{r) and its derivative to be continuous at r = Rci- Instead of /? we will 
often use 5 — {Rci — Rc)/Rc, which characterizes the width of the peak maximum of U{r). Then U{r) is specified 
by two parameters, Ec and 6, with [24| 



Rc 



a{2 + 3S) 



1 



2Ec{l + 5Y' 
Rci = i?c (1 + 5), Eci 



5{2 + 36)' 



U{Rci) = E, 



c 



2 + 25 
2 + i5' 



(4) 



The potential U{r) passes through zero at r = Rco — Rc{^ — P^^^^)] its behavior at smaller r is unimportant (in 
our approximation). Realistic models should correspond to /3 ^ 1 (the small-r slope of U{r) should be sharp; Rco 
should be positive) which translates into (5 <C i (because (3 = 1 corresponds to S = i). 
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B. Sub-barrier energies 



With the potential ^ at E < Ec one has 

=5o exp*(S), (5) 

where is taken analytically and has a rather complicated form given by Eq. (9) in Ref. [l^l- Notice that in our 

semi-classical approximation 'i'(Ec) — iirrjc, where rjc — ri{Ec) is the Sommerfeld parameter at i? = Ec- 

In this paper we propose a simplified expression for S{E) at E < Ec- Let us recall that at i? < Ec is 

accurately approximated [24] by the Taylor expansion '^{E) = go + giE + g2E^ + - . - The explicit expressions for the 
expansion coefficients t^o, ffi, and 52 in terms of Ec and 6 are given by Eqs. (14) and (15) of Ref. [IJ]. At i? < Ec we 
suggest the approximation 

•^{E)=go+giE + g2E^, (6) 
32 = (2^77c - go - giEc)lEl. (7) 

The expressions for go and gi will be taken from Ref. [2^] , while 52 is introduced in such a way to satisfy the condition 
'^{Ec) — "^.nrjc- Our polynomial approximation ^ is much simpler than the exact expression (9) of Ref. [IJ] but 
gives nearly the same accuracy in approximating all S'(£')-factors considered in Sec. |lll Nevertheless, it is still not 
very convenient in applications because Eqs. (14) and (15) of Ref. [13] for go and gi are rather complicated. Therefore 
we further simplify these equations by adopting the limit oi S 1, which is sufficient for applications in Sec. |IVl In 
this limit we have t2JI 
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§^ (8-7rV2J-25), (8) 
Ec 




.9i--Wi?rU-'^v25-,5 . (9) 



Thus, Eqs. (H])-® give simple and practical expressions for S{E) at E < Ec- In our model, S{E) is determined by 
the three parameters, Ec, S, and So', Ec and S determine the shape of the potential U{r); So specifies the efficiency 
of fusion reaction. 

C. Contribution of £ > waves 

The proposed model is phenomenological, being based on s-wave semi-classical tunneling through a spherical poten- 
tial barrier U (r). The actual reaction cross section a{E) contains contribution of different £ waves, with £ = 0, 1, 2, . . . 
For a given multipolarity i, we have a quantum mechanical scattering problem of two nuclei moving in an effective 
potential 



Kff(r) = U{r) + Vi{r), V,{r) = (10) 



where Vg{r) is the centrifugal potential. Then the cross section is (e.g., Refs. [29l. Isoj) 



o{E) = ^Y.^2i^l)UE)P,{E), (11) 

where k is wave-number {E = fi^k'^ /2fi), Ti{E) is the transmission coefficient, and Pe{E) is the fusion probability for 
the penetrating wave. At the low energies of our interest, one traditionally assumes Pi{E) ~ 1. 

Let us return to our S'(£')-model at i5 < Ec- In accordance with pT|) . we have S{E) = J2e S'^^^E), where S'^^^E) 
is an ^-wave contribution to S{E) to be calculated using the effective potential Ves{r). At E < Ec, all i waves refer 
to subbarrier motion, and the transmission coefficient Ti{E) decreases evidently with the growth of ^. Then S{E) can 
be written as 

S{E)^S^'\E)J{E), J{E) = l + Y,{2t+l)^^ (12) 
Here, S^^\E) is the s-wave contribution, and the sum in J{E) is the correction due to waves with t>l. 
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A simple analysis shows that S{E) is typically determined by several lowest i a.t E < Ec for reactions considered 
in Sec. [Ill The correction factor J{E) appears to be essentially higher than 1 (values ^ > are important) but is a 
slowly varying function of E. In other words, the transmission coefficients Ti{E) at these £ are similar functions of 
E. Their main energy dependence is the same as for s-wave, S^°''{E). A crude estimate at E noticeably below Ec 
gives J{E) ^ 1 + y/Ec/Eo, where Eq — fi,^/(2/xi?^) is the characteristic quantum of centrifugal energy (typically, 
Eq ^ Ec)- To simplify the model, a.t E < Ec we suggest the approximation 

S{E) = SosJo exp(go + 9iE + 52^'), (13) 
Jo = l+joVEc/Eo. (14) 

Here, Sqs is the s-wave contribution to Sq, J{E) is approximated by energy-independent constant Jq (that is specified 
by a constant jo); Jq (or jo) can be treated as a parameter which characterizes the importance of higher multipolarities 
£ > 0. One often states in quantum mechanics that the main contribution to scattering at low E comes from £ — 0. 
This statement does not apply to the potential considered here, which has strong attraction at r = 0. In this case, 
higher £ are important even at very low E. 

Finally, at E = Ec, the s-wave barrier is removed and we have Tg{Ec) = 1 in Eq. pip . Then the s-wave partial 
cross section is (Jo{Ec) = 7r/fc^, where kc is the wave- number referring to £^ = Ec- From this (To{Ec) and Eq. ^ 
we immediately get: 

So, = — = 0.6566 ^f^^ MeV barn. (15) 
2^ A1A2 

Therefore, Sqs is specified by Ai and A2. The astrophysical 5- factor at _E < Ec is determined by the parameters Ec 
and P (or 6) of the U{r) potential and by the factor Jq (or jo). 



D. Above-barrier energies 



At £^ > Ec the effective barrier is transparent for some \ow-£ waves. In this case we adopt the simplest barrier- 
penetration model with T({E) = 1 if the Vcs{r) barrier is transparent at a given E, and with Ti{E) — otherwise. 
The cross section is then given by Eq. PT|) with Ti = 1 and Pi = 1, where the sum is taken from ^ = to some 
maximum £oiE) at which Vcs{r) becomes classically forbidden. To be transparent at lower £, the Vcs{r) potential 
should have a pocket (with a local minimum) and a barrier (with a local peak) at r < Rc; let r = vq < Rc be the 
peak point. It is well known that in this case a.t E > Ec the cross section becomes (e.g., Ref. ,29j ) 

'y{E)^^{£o{E) + l)\ (16) 

With increasing E, the range of ^ < £q{E) widens. In the spirit of semi-classical approximation, at £^ > Ec we can 
treat £o{E) as a continuous variable. As long as E is not too much higher than Ec, the local peak point rg of T4ff(r) 
is close to Rc- To locate this peak, it is convenient to introduce x — {Rc — r)/Rc instead of r and linearize Vi{r) in 
terms of x, keeping constant and linear terms. The peak occurs at x = xq ^ (£ + l)£Eo/{EcP) at which dVcs/dx — 0. 
The peak height is Veff (ro) = Ec + {£ + 1)£Eq + {£ + l)^£'^E^/{Ecl3). At a given E the barrier becomes classically 
forbidden when Vce{fo) — E, which gives a quadratic equation for {£0 + l)£o- Solving this equation, we have 



(fo + wo^^^Sl!±MZM^^ (17) 

Now we can easily find £q{E). Substituting it into Eq. (fT6|) gives a closed expression for <t{E) at E > Ec- 

However, such a solution seems overcomplicated. It can be further simplified if we notice that at not too high E 
we typically have E — Ec <C /3Ec and expand the expression containing square root. This gives 

the second term in the parentheses being a small correction. 

One has £o{E) S> 1 very soon after E exceeds the barrier Ec; then {£0 + l)£o ~ £l- On the other hand, at E = Ec 
in the adopted approximation we should have £q(Ec) = and ao{Ec) = T^/kc- Naturally, our approximation of 
continuous £q{E) is inaccurate just near the threshold (E — )■ Ec), where s-wave nuclear collisions proceed above the 
threshold, while higher-^ collisions operate either in the subbarrier regime or in the transition regime (emerging from 
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under respective barriers with growing E). It is a complicated task to describe accurately (j{E) and S{E) at E Ec- 
In order to preserve the simphcity of our S{E) model we propose writing the cross section a{E) at E > Ec as 

<j{E) = ^ ^y^{E) + J^, (19) 

where y{E) is given by Eq. (IT5|) and Jo by Eq. ([Ti)) . Then at i? = Eq this equation matches the cross section model 
proposed for subbarrier energies (Sec. IIIIBp while at i? — Ec JqEq it reproduces the semi-classical cross section 
(HH). 

Translating the cross section ([T9l) into S{E), at E > Ec we finally have 



SiE) = Sos exp(2^77) ^y^E) + Jl (20) 

At E ~ Ec ^ JqEq the leading term in the astrophysical 5-factor is S{E) « Sqs exp(27r?7) (£■ — Ec)/Eq. In the 
previous work we described 5(i?) by a phenomenological expression (Eq. (6) in Ref. [2J|). It reproduces this leading 
term at JqEq <C £' — Ec Ec (see the term containing in Eq. (6) of Ref. [l^l) but diverges from it at higher E. 
Although we do not intend to propose an accurate approximation of S{E) at high energies (few barrier energies Ec 
and higher), we remark that the phenomenological expression (6) for S{E) at E > Ec in Ref. 24] is now replaced by 
Eq. (l20l) which reproduces the correct semi-classical behavior of S{E) above the barrier. Our new expression ((20)) is 
self-sufficient, it contains no extra input parameters in addition to those introduced for describing S{E) at E < Ec 
in Sec. ImB] 



E. How to calculate S{E) 

For any fusion reaction, S{E) is determined by three input parameters, Ec, P (or 5, see Eq. ^) and Jq (or jg, see 
Eq. ([Ml)). The first two parameters specify the shape of the effective potential U{r), Eq. (|31); Jq (or jo) takes into 
account the contribution of higher i = 1,2,... to S{E). At E < Ec the astrophysical factor S{E) is given by Eq. 
([T^. and at E > Ec it is given by (|20p . The factor 5*0^ entering these equations is defined by (fTS)) . The parameters 
go, gi, and g2 in Eq. ([13]) are given by ([5]), (O, and ([7]). 

By construction, our S{E) model can be accurate at i? up to a few Ec- Although the model is generally based 
on first principles, it is phenomenological in a narrow energy range near the barrier {\E — Ec\ ^ JqEq). Our current 
S{E) model is simpler than its previous version [^J] and contains 3 input parameters instead of 4 in [2J]. We will see 
that the current model is more accurate than the previous one. 

Let us stress that our new version implies (5^1 (actually, S < 0.1) meaning that the maximum of the effective 
potential U{r) [see Eq. ([3])] is rather sharp. For broader maximum (0.1 ^ S < 1/3), it would be more appropriate 
to use the same Eq. (1^ at i? > Ec but replace our Eq. (|T31) at i? < Ec by Eqs. (5) and (9) of Ref. [2J, setting 
^o = SqsJq in (9). Although this would complicate the expression for S{E) at E < Ec, the fit parameters would 
preserve their meaning. 



IV. FITS 

Let us approximate all calculated S{E) using the analytic model of Sec. IIIII In Ref. we approximated S{E) 
for the 946 reactions involving C, O, Ne, and Mg isotopes with the first version of the model. That approximation 
is sufficiently accurate, but we approximate those data again together with the new data using the elaborated S{E)- 
model to obtain a uniform database with minimum number of input parameters. 

We consider reactions of each type (each line in Table |I|) separately, and apply the analytic model (fT5)) and (I^Hl) to 
every reaction. In this manner we determine 3 fit parameters, Ec, S and jo, for every reaction. For instance, we have 
3 X 120 — 360 parameters for Si+Si reactions. However, we notice that we can set 6 and jo constant for all reactions 
of a given type (for example, S = 0.0409 and jo = 2.8162 for all Si-fSi reactions) without greatly increasing the fit 
errors. Such constant 6 and jo are given in Table HIl Notice that while jo is constant for a given reaction type, Jq is 
given by the scaling relation (ITi|) and differs from one reaction to another. 

Still, we need to specify the barrier height Ec for every reaction. Collectin g th e values of Ec for all reactions of 
each type, we were able to fit them by the same analytic expression as in Ref. [24|: 



Ec = a/Ri2, Ri2^R + Ai?i 1^1 - Aiol -I- Ai?2 1^2 - A20I, 



(21) 
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where Aiq = 2Z\ and Aiq = IZ^ are mass numbers of most stable isotopes; Ai?i = Ai?ia at A\ > Aiq\ ARi = ARn, 
at Ai < Aiq; AR2 = A/?2a at A2 > A20; AR2 = Ai?2b at A2 < ^20- This gives 5 new fit parameters R, ARia, Ai?2a, 
ARib, Ai?2f, (also given in Table ITU) for each reaction type, and, hence, 7 parameters in total. 

As in Ref. [24], our fit procedure is based on standard relative deviations of calculated (calc) and fitted (fit) S{E)- 
factors. The absolute value of such a deviation in a point E is ri{E) = |1 — Sfit{E) / Sca.\c{E)\. Fitting has been done by 
minimizing root-mean-square (rms) deviation rj^-ms over all energy grid points for all reactions involved in a fit. Column 
9 of Table HIl lists r^rms for all reactions of a given type over all energy grid points (e.g., over 120 x 379 = 45, 480 points 
for the Si-|-Si reactions). Column 10 presents the maximum absolute value of the relative standard deviations, r^max, 
over all these reactions and points. Root-mean-square values r^rms are reasonably small; they vary from rj^ms ~ 0.07 for 
B-l-B, B-l-C, and B-l-N reactions to iirms ~ 0.28 for Si-|-Si reactions. Maximum relative deviations are larger, reaching 
'/max ~ 1.12 for Si-|-Si. Such large values of 77,„ax hide the proper maximum difference of Sfit{E) to S'caic(-E') in a fit. 

To visualize this difference, in column 11 we list the maximum value of the parameter rj(E) that we define as 
7j{E) = Sm{E)/S,^i,{E) - 1 for SfitiE) > 5caic(-B) and as ?i{E) = Sc.UE)/Sm{E) - 1 for Sfit{E) < S^ME)- Thus 
defined, we have ij = r] a.t Sat > 5'caic but 77 > 77 at 5fit < 5caic. For 77 ^ 1 we always have ij rj, but for 77 > 1 the value 
of ^ can be much larger than ij. One can see that Tymax + l is the maximum value among ratios Sm/Scsdc (at S^t > S'caic) 
and Scale/ Sfit (at S^t < S'caic) in a fit sample. For our Si-|-Si reactions we have ?7max + 1 ~ 5.40. This relatively large 
difference between S'caic and Sm occurs at one energy point in one of the 120 Si-|-Si reactions. Specifically, it happens 
for the ■**Si-F**Si reaction (at £; = 12.1 MeV) which involves very neutron-rich nuclei. Although Tymax + l is quite large, 
our fits are still acceptable because the expected nuclear physics uncertainties associated with calculating S{E) for 
such nuclei are larger (SecHT]). Notice that r?max and f^max in Table|TT]can occur for the same reactions and in the same 
energy points, and they can even coincide there. This happens, for instance, for B-l-N reactions (?7max = ??max = 0.37 
for -'^^B-I-^^N a.t E = 6.1 MeV). However, they can also occur for different reactions and at different E (for instance, 
in case of Si-|-Si we have rymax = 1-12 for the ''°Si-|-''°Si reaction a.t E = 14.2 MeV). We could have changed the fit 
algorithm and obtain somewhat smaller r^max (for instance, by minimizing rjj-^g instead of Tyrms) but it would give 
somewhat larger ?7rms- 




E [MeV] 

FIG. 1. (Color online) S-factors for six C-|-Si reactions Ai + A2 {A\ and A2 are given in parentheses). Filled dots are original 
calculations (on a rarefied grid of E points); solid lines are our fits (Table lll)l : open dots show the fit values of Ec- 

A graphical comparison of the previous fits with calculations for many reactions is presented in Figs. 3 and 4 of 
Ref. ^^J]. The comparison with the new fits looks the same. For instance, in Fig. [T] we compare the fitted S{E) 
(solid lines) with calculations (dots on a rarefied grid of iJ-points) for six C-|-Si reactions. The reactions are labeled 
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by (^1,^2), where Ai and A2 are mass numbers of C and Si isotopes, respectively. The lower line corresponds to 
the "'^"C+^^'Si reaction involving the lightest (proton rich) nuclei from our collection (Table HI). The upper line is for 
the ^'*C+''^^Si reaction involving our most massive (neutron rich) C and Si isotopes. Other lines refer to intermediate 
cases, including the -'^^C+^^Si reaction of most stable nuclei. We see that the ^-factors for C+Si reactions vary over 
many orders of magnitude but the fit accuracy remains acceptable. 

The fit errors for the new S'(i?)-model are generally lower than for the previous one [24j. For instance, for the 0+0 
reactions we now have the maximum relative fit error Tymax ~ 0.47 and the rms relative fit error 771ms ~ 0.13 while in 
Ref. |24| we had r^max ~ 0.50 and ryims ~ 0.14 (and we now have 7 fit parameters instead of 9). The fit accuracy is 
especially improved for reactions involving lower- Z nuclei. 

For some reaction types (e.g., for C+C), the fit errors of individual fits (made separately for every reaction) are 
noticeably lower than for all reactions of a given type. This indicates that the fit formula ((2T|) for Ec can be further 
improved. Also, we expect that the fit quality can be improved by introducing some slowly varying function J{E) 
instead of constant Jq at subbarrier energies in Eq. I|14p. and by elaborating the description of S{E) at near-barrier 
energies \E — Ec\ ^ JqEq. However, we think that the accuracy of the new fits is quite consistent with the quality 
of the present data (Sec.|ll|. The fits give a compact and uniform description of calculated S{E) for many reactions. 
They give reliable S{E) in a wide range of energies E because they are based on first principles. In addition, they are 
convenient for including into computer codes. We warn the readers that the calculated and fitted S{E) do not take 
into account resonances. Therefore, one should add the resonance contribution in modeling nuclear burning which 
involves essentially resonant reactions. 

Let us remark that our analytic model for S{E) can also be used to reconstruct the interaction potential U{r) by 
fitting the S(E) data available from experiment or from calculations. Some examples have been presented in Ref. 
(24| . Another example will be given below. 



V. STUDYING UNCERTAINTIES OF S{E) 

This section illustrates another advantage of our analytic S{E) model - its ability to study possible uncertainties 
of astrophysical reaction rates. 

Clearly, any calculation of S{E) contains some uncertainties. First of all, they can be associated with specific 
theoretical model. In our case (Sec. |ll| they are due to using the Sao Paulo potential, the barrier penetration model 
and the NL3 parametrization for deriving nuclear densities in the RHB theory. These uncertainties influence an 
effective potential U{r) and, hence, S{E). Fitting any given S{E) with our model, one can estimate the effective 
potential U{r) [find /3 and Ec in Eq. ([3])]. Assuming reasonable uncertainties of /3 and Ec and using our S'(£')-model 
again, one can easily estimate the expected range of S{E) variations. 

The 5(i?)-model can also be useful to estimate possible effects of dense matter. Consider, for instance, the heating 
of the inner crust of accreting neutron stars in X-ray transients (the so-called deep crustal heating [ij, [1^ Il3 ) • These 
transients are compact binary systems containing a neutron star and a low-mass companion. The deep crustal heating 
is thought to occur mainly due to pycnonuclear reactions in accreted matter when it sinks in the inner crust under the 
weight of newly accreted material. The heating can power llSjl thermal surface emission of these neutron stars that 
is observed in quiescent states of transients (see, e.g., Refs. OllSl)- Pycnonuclear reactions occur due to zero-point 
vibrations of atomic nuclei in a crystalline lattice. Their physics is described, for instance, by Salpeter and Van Horn 
and applied in later work (e.g. @, Hj25, 31] and references therein). Pycnonuclear reactions occur at high densities 
and involve very neutron-rich nuclei [la . |17| . immersed in a sea of free neutrons available in the inner crust. For 
example, consider the powerful •^^Ne-|-'^'*Ne— )-^*Ca reaction at p « 1.7 x 10^^ g cm~^ in the scenario of Haensel and 
Zdunik [1^]. The fractional number of free neutrons among all nucleons in the burning layer is 0.39. The calculations 
of S'-factors are performed for fusion of nuclei unaffected by dense fluid of free neutrons. However, such a fluid can 
compress the nuclei, modify their interaction potential U{r), and hence S{E). 

Accurate calculations of S{E) under these conditions have not been performed; they are complicated and do deserve 
a special study. Nevertheless, the 5(i?)-model allows us to estimate the range of expected S{E) changes. This is 
illustrated in Fig. [5] which shows the model effective potentials (left panel (a)) and respective astrophysical factors 
S{E) for the '^'^Ne-f-'^^Ne reaction (right panel (b)) without and with possible effects of dense matter. 

The filled dots in the right panel of Fig. [2] are our calculated S{E) on a rarefied grid of energies E and the thick 
solid line is our fit (just like in Fig. [T] for some C-|-Si reactions). The dot-dashed line in the left panel is the effective 
potential U{r) used in calculations. According to the Sao Paulo model, the theoretical effective potential slightly 
depends on E. In Fig. [2] it is taken at = 10 MeV. The thick solid line is our model effective potential U{r), which 
is given by Eq. © and plotted for the parameters Ec = 12.137 MeV and 6 = 0.0441 inferred from the fit (Table|lTl. 
We see that fitting the available (here - calculated) S{E) data with our analytic model allows us to reconstruct U{r) 
with sufficiently good precision. More examples of successful reconstructions are given in Ref. (23 |. The three thin 
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FIG. 2. (Color online) Effective potentials U{r) (left panel (a)) and respective astrophysical S'(i5)-factors (right panel (b)) for 
the ^■*Ne+^*Ne reaction. The dot-dashed line on the left panel is the original Sao Paulo effective potential used in calculations 
(it slightly depends on E and is taken at i? = 10 MeV); dots on the right panel are original calculated S{E) on a rarefied 
grid of S-points. Thick solid lines show fitted S{E) and U{r) reconstructed from these fits. The three thin solid lines present 
appropriate effective potential Vcfi{r) including the centrifugal term, Eq. (|10|) . for £=5, 10, and 15 (from bottom to top). 
Long-dashed, short-dashed, and dotted lines refer to three models 1, 2, and 3, which illustrate possible effects of dense matter 
on the '^'*Ne-|-^''Ne reaction (see text for details). 



solid lines in Fig. [5] illustrate our discussion on the contribution of ^ > waves to nuclear fusion (Sec. IIII C|) . These 
solid lines represent the reconstructed effective potential, T4ff(r), including the centrifugal terms £ —5, 10, 15 (bottom 
to top). A few lowest (-waves penetrate the barrier almost with the same efficiency as the s-wave and contribute to 
S{E) even at low E. 

Other curves in Fig. [5] demonstrate possible effects of dense matter. Nuclear reaction rates are proportional to the 
value of S{E) at a typical reaction energy E. In pycnonuclear reactions the energy E is rather low [B*], so that the 
reaction rate is actually determined by 5(0). It is reasonable to expect that the presence of free neutrons between the 
reacting nuclei broadens and/or lowers the maximum of U{r). In Fig. [2] we consider three models of this phenomenon 
labeled as 1, 2, and 3. Model 1 (long-dashed lines) assumes extra broadening of the U{r) peak at the same height 
[Ec = 12.137 MeV as before but larger 6 = 0.1). The barrier U{r) becomes thicker (left panel) which lowers S{0) by 
about eight orders of magnitude. The reaction rate will be strongly suppressed which may cause [l^l a delayed ^"^Ne 
burning after accretion stops in an X-ray transient. Model 2 (short dashed lines) assumes the same curvature of the 
U{r) peak {6 — 0.0441) as for the initial potential but lower maximum, Ec = 11 MeV. The lower barrier is naturally 
more transparent, which increases the initial value of 5(0) by about three orders of magnitude; the '^'*Ne burning will 
react quicker to variations of accretion rate. This possibility has also been considered in Ref. [32| . Finally, model 3 
(dotted lines) assumes that the medium effects simultaneously lower and broaden the barrier {Ec = 11 MeV, 5 = 0.1). 
The lowering makes the barrier more transparent while the broadening makes it more opaque. In this example, the 
broadening wins so that S{0) is about two orders of magnitude smaller than the initial value. Therefore, we may 
expect that the medium effects can greatly enhance or suppress pycnonuclear reaction rates and the net effect is 
not clear. Also, the presence of free neutrons between the reacting nuclei may change U (r) in such a way that the 
approximation ([3]) becomes poor. In addition, the theoretical expression for the pycnonuclear reaction rate through 
S{E) contains serious uncertainties @, H^l which complicate the problem. Further studies are required to clarify these 
points. 

Let us add that nuclear reaction rates in dense stellar matter (especially, in the cores of white dwarfs and envelopes 
of neutron stars [1|) can be greatly affected not only by the transition to pycnonuclear burning regime but also by 
plasma scree ning of the Coulomb interaction. The plasma effects were described by Salpeter and Van Horn [BJ (also 
see [gI. [tI. [25l. 1311133} and references therein). They modify the interaction potential U{r) but mainly at sufficiently 
large r, typically larger than nuclear scales, whereas we discuss the nuclear physics effects which influence U (r) at 
smaller r. It is commonly thought that the plasma physics and nuclear physics effects are distinctly different and can 
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be considered separately. However, as we noticed in [24|, in dense and not very hot stellar matter both effects become 
interrelated and should be studied together. 



VI. CONCLUSIONS 

We have performed new calculations and created a large database of the 5- factors for about 5,000 fusion reactions 
(Sec. HI]) involving various isotopes of Be, B, C, N, O, F, Ne, Na, Mg, and Si located between proton and neutron drip 
lines. Note that the drip lines are obtained in the RHB calculations with the NL3 parametrization. The S-factors were 
calculated using the Sao Paulo method and the barrier penetration model with nuclear densities obtained in the RHB 
calculations. We have elaborated and simplified (Sec. IHII) our model 0| for describing the astrophysical S'-factor as 
a function of center-of-mass energy E of reacting nuclei for non-resonant fusion reactions. Our main results are: 

• For any reaction, we present S{E) in a simple analytic form in terms of three parameters (Sec. IHI Ep . They 
are Ec, the height of the Coulomb barrier; 6 that describes the broadening of the peak of the effective barrier 
potential U{r); and jo that measures the contribution of > waves at subbarrier energies. Analytic fits are 
expected to be sufficiently accurate for energies below and above Ec (up to a few Ec)- 

• We succeeded to fit all our S{E)-da.ta. with only 7 fit parameters for any group of reactions involving isotopes 
of the same elements fSec. llVi Table HIl) . The fit accuracy is well within estimated nuclear-physics uncertainties 
of calculated S{E) (Sec. [HI- 

• In this way we obtain a simple, accurate, uniform and ultracompact database for calculating S{E); the instruc- 
tions for users are given in Sec. IHI El It is easy to implement the database into computer codes (especially in 
network-type ones) which calculate nuclear reaction rates and simulate various nuclear burning phenomena in 
astrophysical environment. 

• In comparison with our previous S{E) model [23|, the present version is simpler, more accurate and reliable. 
We have simplified the analytic expression for S{E) at subbarrier energies (Sec. IIIIB[) . We have clarified the 
contribution of ^ > waves to S{E) and introduced the parameter ( Jq or jo) that accounts for this contribution 
(Sec. IHI C|) . We have replaced a phenomenological analytic expression for S{E) at energies above the barrier by 
a rigorous semi-classical expression (Sec. IHI Dl) . These modifications have allowed us to reduce the number of 
fit parameters (now 3 parameters instead of 4 for any given S{E), and 7 parameters instead of 9 for any group 
of reactions involving isotopes of the same elements) . 

• We have discussed (Sec. |V| the possibility of using our model for estimating uncertainties of S{E) values and 
for studying the effects of dense matter on S{E). For illustration, we have analyzed the range of variations of 
S{E) due to in-medium deformations of interaction potential U{r) for the pycnonuclear '^^Ne-|-^''Ne reaction in 
the inner crust of an accreting neutron star (with the conclusion that the variations can reach several orders of 
magnitude) . 

As detailed in Ref. [l^], the analytic S{E) model is practical for describing large uniform sets of S{E) data. The 
parameters of the model can be interpolated from one reaction to another which can be useful in the case when 
some S{E) data are not available. The functional form of our analytical S{E) is flexible enough to describe different 
behaviors of S{E) at low energies [l^]. Fitting a given S{E) (computed or measured in laboratory) with our analytic 
model can be used to reconstruct the effective potential U{r) (Sec. [V) also see Ref. [l^l). There is no doubt that the 
present model can be improved (Sec. IHI El Ref. [23]), particularly, by complicating the interaction potential U{r). 
However, this would require the introduction of new parameters which would complicate the model. This will improve 
the description of S{E) in some particular cases but the most attractive features of the model - its simplicity and 
universality - would be lost. 
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TABLE I. Fusion reactions (Ai, Zi) + {A2, Z2) under consideration. See text for details 



Reaction 


Ai 


A2 


F, 


Number 


type 






MeV 


of cases 


Be+Be 


8-14 


8-14 


15.9 


10 


Be+B 


8-14 


9-21 


16.9 


28 


Be+C 


8-14 


10 24 


16.9 


32 


Be+N 


8- 14 


11-27 


17.9 


36 


Be+O 


8-14 


12-28 


18.9 


36 


Be+F 


8-14 


17-29 


18.9 


28 


Be+Ne 


8-14 


18-40 


19.9 


48 


Be+Na 


8-14 


19-43 


21.9 


52 


Be+Mg 


8-14 


20-46 


22.9 


56 


Be+Si 


8-14 


24-52 


23.9 


60 


B+B 


9-21 


9-21 


15.9 


28 


B+C 


9-21 


10-24 


16.8 


56 


B+N 


9-21 


11-27 


17.8 


63 


B+O 


9-21 


12-28 


18.8 


63 


B+F 


9-21 


17-29 


18.8 


49 


B+Ne 


9-21 


18-40 


19.8 


84 


B+Na 


9-21 


19-43 


21.8 


91 


B+Mg 


9-21 


20-46 


22.8 


98 


B+Si 


9-21 


24 52 


23.8 


105 


C+C 


10-24 


10 24 


17.9 


36 


C+N 


10-24 


11-27 


19.8 


72 


C+F 


10-24 


17-29 


20.8 


56 


C+0 


10-24 


12-28 


17.9 


72 


C+Ne 


10-24 


18-40 


19.9 


96 


C+Na 


10-24 


19-43 


21.8 


104 


C+Mg 


10 24 


20 46 


19.9 


112 


C+Si 


10 24 


24-52 


24.8 


120 


N+N 


11-27 


11-27 


17.8 


45 


N+O 


11-27 


12-28 


19.8 


81 


N+F 


11-27 


17-29 


20.8 


63 


N+Ne 


11-27 


18-40 


21.8 


108 


N+Na 


11-27 


19-43 


21.9 


117 


N+Mg 


11-27 


20-46 


21.9 


126 


N+Si 


11-27 


24-52 


24.8 


135 


0+0 


12-28 


12-28 


19.9 


45 


O+F 


12-28 


17-29 


21.8 


63 


0+Ne 


12-28 


18-40 


21.9 


108 


0+Na 


12-28 


19-43 


21.8 


117 


O+Mg 


12-28 


18-46 


21.9 


126 


0+Si 


12-28 


24-52 


24.8 


135 


F+F 


17-29 


17-29 


19.8 


28 


F+Ne 


17-29 


18-40 


21.9 


84 


F+Na 


17-29 


19-43 


24.9 


91 


F+Mg 


17-29 


20-46 


24.9 


98 


F+Si 


17-29 


24-52 


29.9 


105 


Ne+Ne 


18-40 


18-40 


21.9 


78 


Ne+Na 


18-40 


19-43 


29.9 


156 


Ne+Mg 


18-40 


20-46 


24.9 


168 


Ne+Si 


18-40 


24-52 


29.8 


180 


Na+Na 


19-43 


19-43 


21.8 


91 


Na+Mg 


19-43 


20-46 


29.9 


182 


Na+Si 


19-43 


24-52 


37.9 


195 


Mg+Mg 


20-46 


20-46 


29.9 


105 


Mg+Si 


20-46 


24-52 


39.8 


210 


Si+Si 


24-52 


24-52 


39.8 


120 
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TABLE II. Fit parameters of S{E) for reactions {Ai, Zi) + {A2, Z2) 



Reaction 
type 


R 
fm 


fm 


Ai?2a 

fm 


AJ^16 
fm 


fm 




5 


io 




^max 




Be+Be 


7.5010 


0. 


2480 


0. 


.2480 


0, 


.1557 





.1557 


0. 


.0330 


0. 


.7453 


0. 


.08 


0, 


.37 


0, 


.59 


Bc+B 


7.5065 


0. 


2547 


0. 


2223 


-1, 


.7469 





.0635 


0, 


.0370 


0, 


.7814 


0. 


08 


0, 


.40 


0, 


.54 


Be+C 


7.6982 


0. 


2543 


0. 


1877 


2, 


.0844 


0, 


.0012 


0, 


.0441 


0, 


.7349 


0. 


09 


0, 


.47 


0, 


.67 


Be+N 


7.9324 


0, 


,2523 


0, 


.1546 


7, 


.7591 


-0 


.0233 


0, 


.0484 


0, 


.7446 


0. 


.10 


0, 


.41 


0, 


.70 


Be+0 


8.0708 


0. 


.2507 


0. 


.1346 


-0, 


.0738 


-0, 


.0271 


0, 


.0509 


0, 


.7699 


0. 


.11 


0, 


.44 


0, 


.80 


Be+F 


8.1585 


0. 


.2510 


0. 


.1201 


-0. 


.6214 


-0, 


.0272 


0. 


.0509 


0. 


.8117 


0. 


.12 


0. 


.48 


0. 


.82 


Be+Ne 


8.1485 


0. 


.2510 


0. 


.1270 


-1 


.5671 





.0007 


0. 


.0513 


0. 


.8697 


0. 


.13 


0. 


.59 


1 


.26 


Be+Na 


8.2139 


0. 


.2494 


0. 


.1161 


-1 


.3477 


-0 


.0059 


0, 


.0505 


0, 


.9065 


0. 


.13 


0, 


.52 


1 


.10 


Bc+Mg 


8.2734 


0, 


.2481 


0, 


.1067 


-0, 


.5794 


-0 


.0128 


0, 


.0499 


0, 


.9528 


0, 


.14 


0, 


.57 


1, 


.30 


Bc+Si 


8.3390 


0, 


.2462 


0, 


.0937 


1 


.8651 





.0126 


0, 


.0481 


1, 


.0538 


0, 


.16 


0, 


.62 


1, 


.26 


B+B 


7.6600 


0, 


.2175 


0. 


.2175 





.0511 





.0511 


0, 


.0459 


0, 


.7772 


0. 


.07 


0, 


.51 


0, 


.51 


B+C 


7.8155 


0. 


.2155 


0. 


.1839 


0, 


.0414 


-0, 


.0148 


0, 


.0511 


0. 


.7919 


0. 


.07 


0. 


.40 


0, 


.43 


B+N 


8.0004 


0. 


.2132 


0. 


.1522 


0, 


.0334 


-0, 


.0359 


0. 


.0523 





.8523 


0. 


.07 


0. 


.37 


0. 


.37 


B+0 


8.1037 


0. 


.2117 


0. 


.1331 





.0338 


-0 


.0378 


0, 


.0520 


0, 


.9001 


0. 


.08 


0, 


.47 


0, 


.47 


B+F 


8.1700 


0. 


.2114 


0. 


.1188 


0, 


.0349 


-0, 


.0379 


0, 


.0502 


0, 


.9575 


0. 


.09 


0, 


.48 


0, 


.51 


B+Ne 


8.1755 


0. 


.2088 


0. 


.1244 





.0113 





.0098 


0, 


.0501 


1, 


.0353 


0. 


.12 


0, 


.57 


1, 


.08 


B+Na 


8.2418 


0, 


.2075 


0, 


.1137 





.0082 





.0160 


0, 


.0492 


1, 


.0885 


0, 


.11 


0, 


.56 


0, 


.85 


B+Mg 


8.3033 


0, 


.2064 


0, 


.1047 





.0061 





.0211 


0, 


.0485 


1, 


.1431 


0, 


.12 


0, 


.58 


1, 


.03 


B+Si 


8.3977 


0. 


.2048 


0. 


.0921 


0, 


.0035 


-0, 


.0201 


0. 


.0472 


1 


.2073 


0. 


.14 


0. 


.59 


1 


.05 


C+C 


7.8843 


0. 


.1836 
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